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Abstract. Suppose a graph-directed iterated function system consists of maps /e 
with upper estimates of the form d(/e(a;), /e(i/)) < r^d(x,y). Then the fractal di- 
mension of the attractor of the IFS is bounded above by the dimension associated 
to the Mauldin— Williams graph with ratios r^- Suppose the maps /e also have lower 
estimates of the form d[fe{x),fe(y)) > r'^ d(x,y) and that the IFS also satisfies the 
strong open set condition. Then the fractal dimension of the attractor Ky of the 
IFS is bounded below by the dimension associated to the Mauldin— Williams graph 
with ratios r^. When re = r^, then the maps are similarities and this reduces to the 
dimension computation of Mauldin & Williams for that case. 



0. Introduction 

Fractal sets may be constructed in many different ways. Barnsley [3] singled out 
the "iterated function system" method: The fractal set K is made up of parts, each 
of which is a shrunken copy of the whole set. Mauldin & Williams [13] provided 
a more general setting, where several sets Ky are involved, each of them is made 
up of parts, and each part is a shrunken copy of one of the parts (the same one or 
a different one). The combinatorics of the way in which the parts fit together is 
described by a directed multigraph. This is described in detail below (Definition 
1.5). In addition to Mauldin and Williams, compare "recurrent iterated function 
system" [3, Ch. X], "Markov self-similar sets" [20], and "mixed self-similar systems" 
[1]; see also [5], [19]. 

In the text [6, §6.4] there is an exposition of the Mauldin-Williams computation 
of the dimension of the attractors for a graph-directed iterated function system 
consisting of similarities. "Similarities" are functions /: 5 — s- T between metric 
spaces that satisfy equations of the form 

(0.1) d{fix)J{y)) =rd{x,y) iova\lx,yeS. 

Here we have used d for the (possibly different) metrics in the two metric spaces. 
This dimension computation involves the spectral radius of some nonnegative ma- 
trices. Next in the text are two "Exercises" [6, (6.4.9) and (6.4.10)] asking for 
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appropriate generalizations when the equahties (0.1) are replaced by inequalities of 
the form 

d{f{x),f{y))<rd{x,y) for all a;, y e 5 

or of the form 

d{f{x),f{y)) > rd{x,y) for all x,y G S. 

In this paper we provide some possible solutions for those exercises. The non-graph 

case (or, in our language, a graph with only one node and several loops) may be 
found in [9, Theorems 9.6, 9.7]. 

We also take this opportunity to make a generalization beyond the setting used in 
[6, §6.4]: wc allow complete metric spaces other than Euclidean space. (Even when 
Euclidean space is the main interest, it is useful to apply the results in ultrametric 
"model" spaces consisting of strings or sequences.) This means that for the lower 
bound estimates wc must use the "strong open set condition" rather than the 
original "open set condition" (this has been discovered since the publication of [6] : 
for example [2], [16], [17]). 

The paper is arranged as follows. Section 1 reviews the notation for graph- 
directed iterated function systems. Section 2 reviews the Perron-Frobenius theory 
to be used. Section 3 contains the proofs. Section 4 describes an example. 

This paper is based in part on a portion of the dissertation [11] of the second 
author, written at The Ohio State University under the direction of the first author. 



1. The Setting 

We will describe here: the "fractal dimensions" that will be used, the "graph- 
directed iterated function systems" and the fractal sets Ku that they define, and 
the "string models" that will be used in the investigation. The notation will 
mostly follow [6], especially Section 4.3 on graph self-similarity. 

BVactal Dimensions. The first dimension we consider is the box-counting dimen- 
sion. See [9] for a more complete discussion. 

Let 5 be a metric space, and let C 5 be a totally bounded set. For each 
S > 0, define Ns{K) to be the smallest number of sets of diameter < S that can 
cover E. Since we have postulated that K is totally bounded, Ns{K) is a finite 
natural number for every S > 0. Unless K is a. finite set, Ns{K) increases to oo as 5 
decreases to 0; the rate at which Ns{K) increases will tell us something about the 
size of the set K. 

(1.1) Definition. The upper box-counting dimension or Bouligand dimen- 
sion of K is 

dimB-fV = limsup — . 

5^0 - log 

The lower box-counting dimension of K is 

dini TiK = liminf — ^— -^-^-r^. 
" s^o - log S 



Next we will review the definition of the Hausdorff dimension. See for example 
[6, p. 149], [9, p. 28]. 
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(1.2) Definition. Lot S* be a metric space, let K C S, and let (5 > be a positive 
number. A collection C of subsets of S is S-Rne iff diam A < S for all A G C. The 
collection C is a cover of K iS K C{J^^^^ A. 

(1.3) Definition. Let 5 be a metric space, let K C S, and let s > 0. For ^ > 0, 
define 

) = inf ^ (diam 

where the infimum is over all countable (5-fine covers C of i^. Define 

:K'{K) = sup ■Kl{K) = \im%l{K). 

6>0 S^O 

The Hausdorff dimension of the set K is 

dimn = sup { s : ?{"(ii:) = 00 } = inf { s : K'{K) = } . 

Another fractal dimension (not defined here) is the packing dimension dimp. 

See [7, §1.2] for a discussion. In [6] the packing measure is mentioned only as an 
afterthought, but in [7] its role is nearly as important as the Hausdorff measure. 

Many other fractal dimensions may be found in the literature. 

It is known ([7, Corollary 1.3.5] and [9, Lemma 3.7]) that all of the dimen- 
sions mentioned here (and many other fractal dimensions, as well) are between the 
Hausdorff dimension and the upper box dimension: 

dimn K < dimK < dimB^sT, 

where "dim" is any of dim^ , dimp , dim ^ , dimR . 

Graphs and Iterated Function Systems. First, a directed multigraph {V,E) 
should be fixed. The elements v G V are the vertices or nodes of the graph; the 
elements e G E are the edges of the graph. For u,v G V, there is a subset Eu^ 

of E, known as the edges from u to v. Each edge belongs to exactly one of these 
subsets. We will sometimes write Eu = [JyEuv, the set of all edges leaving the 
vertex u. 

We will often think of the set E a.s a set of "letters" that label the edges of the 
graph, so we will talk about "words" or "strings" made up of these letters. A path 
in the graph is a finite string a = 6162 • • • e/j of edges, such that the terminal vertex 
of each edge is the initial vertex of the next edge e^+i. The initial vertex of a is 
the initial vertex of the first letter ei and the terminal vertex of a is the terminal 
vertex of the last letter Cfc. Wo write E^ for the set of all paths of length k that 
begin at u and end at v; and Ei^'' for the set of all paths of length k that begin at 
u; and Eu^ for the set of all finite paths of any length that begin at u; and E^*'^ 
for the set of all finite paths. By convention, E^^ consists of a single "empty path" 
A„ of length zero from node u to itself. 

If a e E^ and (3 G , then we write afi for the path made by concatenation 
of the two strings, so that a/3 G Euw ■ 

A partial order may be defined on E'<*^ as follows: write a < /3 iff a is a prefix 
of /?: that is, (3 = a'y for some path 7. With this ordering, each set Eu^ becomes 
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a tree. The root of the tree Eu is the empty path A„. The entire set i?^*^ is a 
disjoint union of trees (a forest). Two strings are incomparable iff neither is a 
prefix of the other. 

A nonempty path that begins and ends at the same node is called a cycle. 

We will assume that the graph {V, E) is strongly connected, that is, there is 
a path from any vertex to any other, along the edges of the graph (taken in the 
proper directions). We will also assume that there are at least two edges leaving 
each node. (In [8] may be found a method to eliminate this assumption.) 

(1.4) Definition. An iterated function system (or IFS) directed by the graph 
(y,E) consists of: one complete metric space Sy for each vertex v V and one 
function fg-.Sy^ Su for each edge e e Euv 

To avoid trivialities, wc will normally assume that each metric space Sv has 
at least two points. For this sort of graph-directed arrangement see [13], [6], and 
references in both. Note that in some places, such as [19], the direction of the edges 
of the graph arc the reverse of the convention used in this paper; this may lead to 
consideration of the transpose of our matrix, but of course that does not affect the 
spectral radius. 

(1.5) Definition. An invariant set list for the iterated function system (/e) 
consists of one nonempty set C for each v &V such that 

for all u gV. 

If the functions satisfy appropriate conditions, then it may be proved that 
there is a unique invariant set list of nonempty compact sets. For example [9, 
Theorem 9.1], [6, Theorem 4.3.5]. 

The Models. We will use some "string models" in our investigation of these in- 
variant sets. Write £"!"''' for the set of all infinite strings, using symbols from E, 
where the initial vertex of the first edge is u and the terminal vertex of each edge is 
the initial vertex of the next edge. These sets are naturally compact metric spaces. 
For each finite string a G E^*\ the cylinder [a] is the set of all infinite strings 
a e E^'^^ that begin with a. Then the set { [a] : a G Eu'^ } is the base for the 

topology on E^\ For a G E^^ and a positive integer fc, the restriction a\k is 
the finite string made up of the first k letters of a. The same notation a \k is used 

for finite strings a when k is less than the length of a. As a special case of this: if 
a has length fc, then the parent of a is a~ = a\{k — 1), obtained by omitting the 
last letter of a. 

Suppose [E. V) directs an iterated fimction system (/e) consisting of contrac- 
tions. Then we may define an "action" of V starting from the maps /g. If 
a = 6162 • • • Gfc e Euv , then the composition = /e^ o /ea o • • • o /cfc is a function 
from Sj) to Su- We will usually write simply q(.t) for fa{x). 

There is a model map /i„ : E^'^ Su for each u, defined so that /i„(ct) is the 
unique element of the set 

00 

f]{a\k)K,,. 

k=l 
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(Of course, here Vk is the terminal vertex for the fcth letter of a.) Then clearly = 
hu\Eif^\, since the sets hu\E'^'^\ satisfy the defining conditions for an invariant 
set list of the IFS (/e). If hu{cr) = x, then we say that the string a is an address 
of the point x. 

Mauldin— Williams Graphs. The computations to be done here can be stated 
in terms of a combinatorial structure that assigns a positive number to each edge 
of the graph. 

(1.6) Definition. A Mauldin— Williams Graph is a directed multigraph {V, E) 
together with a positive number re for each edge e G E. 

If a = 6162 • • - Gfe e E^'^^ is a path, we define r{a) = HiLi'^ei; thus r(a/3) = 
r{a)r{(3) and r(A„) = 1. 

(1.7) Definition. The system (re) is strictly contracting iff re < 1 for all e. 

2. Perron-Frobenius Theory 

The computations will use some information from the Pcrron-Frobcnius theory 
of nonnegative matrices. This material may be found in [10], [15], or [18]. We will 
state here the information that is used below. 

Let A be an n X n matrix. Wo say that A is nonnegative, and write ^ > 0, if 
all entries of A are nonnegative. When A and B are both n x n matrices, we write 
A>BiiA — B>0. Similarly, if x is an n x 1 column vector we write x > if all 
entries arc nonnegative, and we write x > of all entries are positive. 

The spectral radius of an n x n matrix A is the maximum absolute value of 
the complex eigenvalues of A. Write p{A) for the spectral radius of A. 

An nxn matrix A is reducible iff the index set {1, 2, • • • ,n} can be partitioned 
into two nonempty sets /, J such that for all i € I, j € J, the entry of A in row i 
column j is zero. The matrix A is irreducible iff it is not reducible. 

Wc cite here a few results from the Perron-Frobenius theory of nonnegative 
matrices: 

(2.1) If A > is irreducible, then the spectral radius A = p{A) is an eigenvalue of 
A and there is a positive eigenvector x > with Ax = Ax. [10, Theorem 2, p. 53], 
[15, Theorem 4.1, p. 11], [18, Theorem 1.5, p. 20]. 

(2.2) U A> 0, A irreducible, x > 0, x ^ 0, ^x = Ax, then A is the spectral radius 

p{A). [10, Remark 3, p. 63], [15, Theorem 4.4, p. 16], [18, Theorem 1.6, p. 20]. 

(2.3) li A > B > 0, where A is irreducible, then p{A) > p{B). [10, Lemma 2, 
p. 57], [15, Corollary 2.2, p. 38], [18, Exercise 1.9, p. 25]. 

Dimension of a Graph. The "dimension" for a Mauldin-Williams graph is com- 
puted as follows (see [13], [6]). 

(2.4) Definition. Let (V,E, {re)eeE) be a strictly contracting Mauldin-Williams 
graph. For each positive real number s, let M[s) be the matrix (with rows and 
columns indexed by V) with entry 

in row u column v. Let '^{s) = p(^M{s)) be the spectral radius of M{s). Then $ 
is continuous, strictly decreasing, $(0) > 1, lims_>oo ^(s) = 0. The unique solution 
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Si > of <J>(.si) = 1 is the dimension associated to tlie Mauldin- Williams graph 

Note that, as a consequence of the rule for matrix multiplication, if fc € N, then 
the fcth power M{s)^ has entry 

in row u column v. Because of our conventions about empty strings, this is true 
even for A; = 0. 

Suppose that {V,E) is a strongly connected graph. Since M(si) has non-zero 
entry in row u column v whenever there is an edge in Euv, it follows that the matrix 
is irreducible. But 1 is an eigenvalue, so by (2.1) there is a positive eigenvector: 

(2.5) A„ > 0, A„ = E E Y^^u = l- 

veV eeE^^. u&V 

(For convenience, we will write Amin = min^ At, and Amax = max^A^. Thus < 

Amin ^ Amax ^ l*) 

This summation property shows that if we define 

(2.6) iiu{[o]) = r{ay' A^ for all a e E^^J , 
then we have the consistency property 

when a G EuJ- Thus the /i„ extend to Borel measures on the spaces Eu ^ • 

3. The Proofs 

Let iy, E) be a strongly connected directed multigraph, let (S'„)„gy be a family of 

nonempty complete metric spaces, and let { fe)ee b be a family of functions /e : Sy — > 
S'u if e G Euv Suppose each of these functions /e satisfies a Lipschitz condition 

(3.1) d{fe{x),fe{y)) <red{x,y) 

iov x,y e Sv,e G E^v- We assume all < 1 (although [6, Exercise 4.3.9] it is 
actually enough to assume that the system of ratios (rg) is "contracting" in the 
sense that all cycles have ratio < 1)- It follows that there is a unique list of 
nonempty compact sets C Sy satisfying the invariance condition 

Ku= [j fe[Ky\ 
eeE 

for all u £V. For future use, we define 

c^min = min { diam : u gV} , rfmax = max { diam Ku : u gV} . 
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Following Definition 2.4, let si be the dimension associated with the Mauldin- 
Williams graph {V,E, (re)eeB)- Write M{s) and $(s) as in that definition. 
We will show (Theorem 3.5) that the upper box dimensions satisfy 

dime Ky < si 

for all V. (See, for example, [9, §9.3] for a proof of this in the non-graph case and 
in Euclidean space.) 

Now suppose there are also lower bounds of the form 

(3.2) d{f,{x),fe{y))>r',d{x,y) 

for x,y € Sv,e G Euv Let S2 be the dimension associated with the Mauldin— 
WilUams graph (V. E. ('r[.),.eE) ■ Write M'[s) for the matrix and $'(s) for the 
spectral radius in Definition 2.4. 

Under the assumption of the "strong open set condition" , we will show (Theorem 
3.14) that the Hausdorff dimensions satisfy 

dimn > S2 

for all V. (See, for example, [9, §9.3] for a proof of this in the non-graph case, in 
Euclidean space, and in the disjoint case.) 

Certainly r'^ < (since each has at least two points, as we have assumed). 
So M'{s) < M{s) and by (2.3) <^'{s) < $(s) so that S2 < si. The two estimates 
together will give us 

S2 < dimKy < si 

where "dim" is any of the fractal dimensions mentioned above. 

Upper Bounds. For the proof of the upper boimd, we need to estimate the counts 
Ns{Ku). Wc will use "cross-cuts" of our forest of finite strings for this. 

(3.3) Definition. A cross-cut is a finite set T C E'-*^ such that, for every infinite 
string a G E^'^^ there is exactly one n with the restriction a\n gT. 

If T is a cross-cut, we will write T„ = T n Ei*^ and T„„ = T n Ei*J for u,vgV. 
Then, according to the definition, for each u G V, the cylinders {[a] : a € T„ } 

constitute a partition of E^^K 

Let iiu be the measures on the spaces Eu^ as in (2.6). If T is a cross-cut we have 

(3.4) A„ = i^^Ei^^) = ^ „{[a]) = ^ E riarX.. 

(3.5) Theorem. Let fe satisfy upper bounds (3.1). Let Si be the dimension for the 
Mauldin-Williams graph (V, E, (re)eeB) • Then the attractors Ky satisfy dim^Ku < 

Si. 

Proof Let 5 > be fixed. Then 

(3.6) T =\^a:u,v GV,aG E^^^ ,r{a) <S< r(a-) } 
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is a cross-cut (called "first time less than (5"). For each a S T, we have i5rmin < 
r{a) < 6. We may estimate the cardinalities of the sets T„: 

A„ = ^ ^ r{ay'\v>{Sr^iny'Xmin#Tu 
vev aeT^v 

and therefore, for every u G V, we have 

'^min 

But of course the set Ku is covered as follows: 

and the covering sets have diameters 

diam aKy < r{a) diam Ky < ^rfmax- 

Therefore we have 



or (write i] = 5(ii„ax) 



si 



Taking logarithms, dividing by the positive number — logr/, and letting — > 0, we 
get 

— logiV^^ 

dimB-fv« = limsup < Si- U 

-logr? 



Lower Bound: Disjoint Case. Now we will consider the lower bound estimate 
on the dimension of the sets Ky that we get from the inequalities (3.2). Note that 
wc continue to assume upper boimds (3.1) as before, so that the fe are continuous, 
the images aK^ are compact, and so on. (Thus < 1 for all e.) 

Now as before, 1 is an eigenvalue of the matrix M'(s2), so by (2.1) there is an 
eigenvector with positive entries: 

(3.7) a;>o, al^^^t^a;, E^» = i- 

veVeeE„^ uev 

(As usual, we will write A^^^^ = min„ A^ and X[nax = max^ X'y.) If fi[^ are defined by 
= '^'{o'Y'K for all a e Ei*J, then they extend to Borel measures on ^ . 

(3.8) Definition. We will say that the graph-directed iterated function system 
(fe) with attracting sets (Ky) falls in the disjoint case if, for any u,v,v' e V, 
e e Euy, e' e Euy' ,e^e', we have /e [Ky] n /«' [Ky-] = 0. 
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(3.9) Theorem. Let fe satisfy lower bounds (3.2). Let S2 he the dimension for 
the Mauldin-Williams graph {y,E,{r'^)e^E) ■ Suppose the IFS falls in the disjoint 
case. Then the Hausdorff dimensions satisfy diniH Ky, > S2 ■ 

Proof. Because of the disjointness, and because of the compactness of the sets K^, 

there is ?7 > so that dist (/e[_ft'„], /e' [-K^-u']) > for u,v,v',e,e' as in the 
definition. We claim that for any incomparable a, a 

dist [aKy,a'Ky/^ > r'{a~)r}. 

To see this, let 7 be the longest common prefix of a, a' . Because they are incompa- 
rable, they are both strictly longer than 7. (Thus 7 < a~, so r'(7) > r'{a~).) 
There are distinct letters e, e' so that a = je^, a' = 76 Say /3 e E^^l 
/3' e eI^)^,; these strings may be empty. Then l3Ky C and /S'Ky/ C K.u,r, 
since e 7^ e' we have dist{e(3Ky, e' (3' Ky/) > dist{eKyj, e'Kyj>) > rj; so finally 
dist(7e/?if„,7e'/3'i^^') > r'(7)r? > /(a")//. 

Let B C Su he a Borel set. Let S = diam B/r] and let T be the cross-cut "first 
time less than 5" 

(3.10) T = i^a:u€V,v €V,a€ E^^J,r'{a) <S< /(a") } . 

The elements a are all incomparable, and diam B = rjS < rjr'la^), so at most 
one set aKy with a e Tu can meet B. This means h~^[B] C [a] for some a £ T, 

and thus MU/i-M^]) < mUH) = r'iay^K < ^'m..V-''{die.^ SfK 
Now if 6 is a countable cover of Ku by Borel sets, we have 

^ (diam B)- > (^) Y.,.'y{h''[B]) > (^) 

It follows that !K'^ (Ku) > {v'^ /AJnax) > 0- Therefore dimn Ku>S2. □ 

Lower Bound: Strong Open Set Condition. To prove a lower boimd for the 
dimension of an IFS we must of course limit the overlap. This is often done with an 
"open set condition" . In general complete metric spaces (not necessarily subsets of 

Euclidean space) we require the "strong open set condition" . For the strong open 
set condition, and its application in this calculation, we follow [2], [16], [17]. 

(3.11) Definition. Let {V,E) be a directed multigraph, let (S't,)t,gv be complete 
metric spaces, let fe'. Sy ^ Su for e G Eyy be a strictly contracting graph-directed 
IFS, and let (if„)„gy be the attractor list. We say that the IFS (/e) satisfies the 
open set condition iS^, for each v gV there is an open set Uy C Sy such that: 

(a) For all u,v G V, e G Eyy, we have fe [Uy'j C J7„. 

(b) For all u, v,v' eV, e G Eyy, e' G e 7^ e', we have /e [(7„] n/e' [[/„'] = 0. 
We say that the IFS (/e) satisfies the strong open set condition iff there exist 
open sets Uy satisfying (a), (b), and 

(c) For all w G F, we have 'Uy<r\Kyi^ 0. 

Unlike the case where each Sy is Euclidean space, it now no longer follows that 
the images fe [Uy] are open sets. Even if fe is a homeomorphism of Uy onto fe [Uy'j , 
it need not follow that fe [Uy] is open. 

Notice the following consequence of the open set condition. (Of course, (3.2) 
implies that fe is injective.) 
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(3.12) Lemma. Let fe satisfy the open set condition with open sets Uy. Suppose 
the maps /e are injective. If a € EuJ , (3 € E^^, are incomparable, then fa \Uv\ H 

fc.'[Uy,] =0. 

Proof. Let 7 be the longest common prefix of a, a'. Because they are incomparable, 
they are both strictly longer than 7. There are distinct letters e, e' so that a = je/3, 
a' = 7e'/3'. Say /3 G E^I, P' G E^*}^,] these strings may be empty. Now by (a), we 
have /3J7^ C and l3'Uy> C Uyj>. Next, e/?C/„ C fe[U^] and e'l3'U^> C f^.[Uyj>] and 
by (b) these are disjoint. Since the maps fe are injective, so are their compositions, 
so we conclude that aUy = 'yepUy and a'Uy' = "ye' (3'Uyi are also disjoint. □ 

From the uniqueness properties of the invariant sets Ky, it is easy to see that 
Ky C Uy. But Ky C Uy nccd not be true. The strong open set condition provides 

a substitute for this. 

(3.13) Lemma. Let be a strongly connected directed multigraph. Suppose 
the IPS (fe) satisfies the strong open set condition with open sets Uy. Por each 
V gV, there is a cycle C € Eyy such that (Ky C Uy. 

Proof. By (c), there is a point x G Ky f] Uy. Let a G E^"^^ be an address of 
X. Then the singleton {x} is the decreasing intersection of compact sets of the 
form {a\n)Ku^. Now Uy is a neighborhood of the point x, so we conclude by 
compactness that there is n so that {a\n)Ku^ C Uy. Because the graph {V,E) 
is strongly connected, there is a path 7 G Ei*Jy. Let ( = (a ["71)7. Then we have 
-/Ky C Ky^ so that (Ky = {a\n)jKy C {a\n)Ky^^ C Uy. □ 

(3.14) Theorem. Let fe satisfy lower bounds (3.2). Let S2 be the dimension for 
the Mauldin-Williams graph (V,E, {r'e)eGE)- Suppose (fe) satisfies the strong open 
set condition. Then the attractors Ky satisfy dimn Ky > S2. 

Proof. Let s < S2. We will show that dimn Ky > s. Now $'(s) > 1, so the spectral 
radius t = p{M'{s)) is > 1. 

Let Uy be open sets for the strong open set condition. For each v gV, choose a 
cycle € Eyy so that CvKv C Uy. Then 

c = minr'(C^)* 

is a positive number. Let n G N be such that cf" > 1. 

Now we will define a new IFS. First we define a new directed multigraph (V, E). 
The vertices of the new graph will be V as before. The edges from vertex u to 
vertex v will be 

E^y = [aCy:aG E^^ } . 

The corresponding maps are fac,^ = faCv- Note that aC,y[Ky] C K^, so the new 
attractors Ky are subsets of the old attractors Ky. 

We claim that that the new IFS falls in the disjoint case. Indecid, let a(y G Eyy 
and a'(v' G E^y' be distinct edges. Then (yKy C Uy and Cy'Ky/ C Uy'. But a and 
a' have the same length and are unequal, so they are not comparable, and thus by 
Lemma 3.12 we have a(yKy n a'^v'K'y C aUy fl a'Uy' = 0. 
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Figure 4.1. A Julia set 

Now the matrix M'{s) for the new IFS has entry 

in row u column v. But this entry is 

which is the corresponding entry in the matrix cM'{s)". Now by (2.3), 

p(M'(s)) > p(cM'(s)") = ct" > 1. 
This shows (by Theorem 3.9) that dimu Ky > s, so also dimu Ky > s. □ 

4. Example: A Julia Set 

We will consider an example showing how the theorem may be applied to esti- 
mate the fractal dimension of a set. This example is the Julia set for the trans- 
formation z"^ ~ 1/2 of the complex plane C. See, for example, [9, Chapter 14] for 
an explanation of Julia sets in simple cases like this. McMuUen [14] has provided 
an algorithm that allows computation of Hausdorff dimension up to any desired 
accuracy of conformal expanding open maps (so having Markov partition) which 
include this example. 

The Julia set we will consider here is shown in Figure 4.1. See [11] for more 
details of the following computation. 

We divide the set into four parts A, B, C, D, using the four quadrants of the 
plane. Then under the two inverse maps zt-y/w -1-1/2 each of these parts is made 
up of the images of two of the parts. A is made up of images of A and B; B is 
made up of images of C and D; C is made up of images of A and B; D is made up 
of images of C and D. Thus our directed graph will have four vertices and eight 
edges. See Figure 4.2. 

For the upper and lower bounds, we will use a simple estimate based on the 
derivatives. 
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Figure 4.2. The graph 

(4.3) Lemma. Suppose ^{z) — — \/2 maps convex set U C bijectively onto 
convex set W , and ^ W . Write f{w) = + for the branch of the inverse 
that maps W onto U . Let 

m = inf { |z| : z e [/} ; M ^ snp {\z\ : z € U } . 

Then, for any wi,W2 G W, 

\W2 -'Wi\< \f{w2) - f{wi)\ < \W2 ~ Wl\. 

IM III I IjYi I I 



Proof. Note that ^'(z) = 2z, so \ip' {z)\ < 2M for aU zeU. Also f'{w) = l/(2/(u))) 
and if w eW then z = f{w) S U, so \f'{w)\ < l/{2m) for aU w eW. 

Now let wi,W2 G W, and write zi = f(wi), Z2 = f{w2). Since U is convex we 
may integrate along the line segment joining zi to Z2'. 



\W2 - Wll = W{Z2) - V{zi)\ = 



ip' {z) dz 



< { sup \(p'{z)\ ] \Z2- Zl \ = 2M\z2 - Zl\ 



Consequently we have \f{w2) — f{wi)\ > (l/{2AI))\w2 — wi\. Similarly, 



f'iw) dw 



< sup \f'{w)\ ] \w2 - Wi \ = — \W2 -Wl\. □ 

\w€W J '2m 



For the convex sets required here, we begin with the smallest circle that surrounds 
our set, a certain parallelogram inside, and segments on the imaginary axis. We get 
four three-sided regions, each containing the portion of the curve in one quadrant. 
See Figure 4.4. 

The parallelogram is chosen so that for each of the four regions bounded by our 
curves, the images under ±y^w + 1/2 are inside the original regions. (See Figure 
4.5.) 
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Figure 4.4. Inner and outer bounds 




Figure 4.5. Images are inside the regions 



Now we are in a position to estimate the constants for the inequalities. The 
outer circle has radius (1 + a/3)/2, so M = (1 + a/3)/2 and our lower bounds are 
all r'(e) = 1/(1 + \/3). Then 



M'{s) 



a a 

a a 

a a 

LO a a 
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Figure 4.6. Eight regions and their images 

where a = (1 + \/3)~^. The spectral radius is 1 when s = S2 = log 2/ log(l + \/3 ) > 
0.689. 

Similarly, for the closest point we get m = \/9 + 3V3/6, and r(e) = l/(2m) « 
0.7962. Computing the spectral radius as before, we get si < 3.042. 
The (rather poor) estimates are 

0.689 <dmiK < 3.042. 

But that need not be the end of our work. Consider the eight small regions in 
Figure 4.5. Their images, in turn, fall as shown in Figure 4.6. This graph-directed 
iterated function system has 8 nodes and 16 edges. We may work with an 8 x 8 
matrix. [Actually, because of symmetry, we may reduce this to a 2 x 2 matrix, but 
such questions of cfRciency are not dealt with in this paper.] Using the upper and 
lower estimates for these regions, we obtain 

0.735 <dimK < 1.758. 

A considerable improvement, but still not very good. 

If we have a computer to help us, wc may continue to get better approximations. 
When we reach a matrix of size 2^^ x 2^^ our estimate is 

1.069 < dim < 1.077. 

The size of the computation is not as bad as the size of the matrix may suggest: 
each row has only two nonzero entries. 

McMuUen's method [14] computes the dimension for this Julia set as 1.07336 
(according to the computer program on his web site). This agreement provides 
evidence that our methods (and our computer program) are correct. 
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Other Julia sets for the quadratic maps z'^ + c are considered in [11]. For c 
outside the Mandelbrot set, and for c in the main cardioid of the Mandelbrot set, the 
process required to carry out the inner and outer approximations was automated. 
It is hoped that details of this will be published elsewhere. 

Note that the functions + c used in this example are "conformal" in the sense 
used in geometry (except at ^; = 0). So this example is a "self-conformal" set. See 
[12] and [19] for more information on self-conformal fractals. 

For a map f : S ^ T oi metric spaces, one might say that / is conformal iff, 
for every xq G S, there is a positive constant c{xo) such that for every e > there 
is (5 > such that 

d{x,xo) 

whenever d{x, xq) < 5. But as far as we know, the only cases that have been consid- 
ered (aside from similarities in general metric spaces) are differentiable manifolds. 
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